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Abstract 

We propose a new scenario of baryogenesis, in which annihilation of axion domain walls 
generates a sizable baryon asymmetry. Successful baryogenesis is possible for a wide range of 
the axion mass and decay constant, m ~ 10® — 10^® GeV and / ~ 10^® — 10^® GeV. Baryonic 
isocurvature perturbations are significantly suppressed in our model, in contrast to various 
spontaneous baryogenesis scenarios in the slow-roll regime. In particular, the axion domain wall 
baryogenesis is consistent with high-scale inflation which generates a large tensor-to-scalar ratio 
within the reach of future GMB B-mode experiments. We also discuss the gravitational waves 
produced by the domain wall annihilation and its implications for the future gravitational wave 
experiments. 
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I. INTRODUCTION 


Axions may be ubiquitous in nature. Indeed, there appear many axions through com- 
pactihcation of the extra dimensions in the string theory [1, 2], Some of them may remain 
relatively light and play an important role in cosmology such as inflation, dark matter 
and dark energy. In this paper we shall present a new scenario of baryogenesis, in which 
axions play a key role. 

The axion exhibits a shift symmetry, 

a a + C, (1) 

where U is a real transformation parameter. While the shift symmetry keeps the axion 
potential flat at the perturbative level, non-perturbative effects break the symmetry to a 
remnant discrete one. 

Let us suppose that one of the non-perturbative effects gives the dominant contribution 
to the axion potential, which is expressed as 

V{a) ~ f (^1 - cos , (2) 

where m is the axion mass and / is the decay constant. Then, the axion potential has a 
series of N (approximately) degenerate vacua, where the precise value of N depends on 
the details of the UV theory.^ If the axion is lighter than the Hubble parameter during 
inflation, it acquires quantum fluctuations which extend beyond the Hubble horizon. For 
sufficiently large quantum fluctuations, some of the N vacua might be populated, which 
results in domain wall formation after inflation. The domain walls are cosmologically 
problematic, and so, they must annihilate before dominating the Universe. This is possible 
if the degeneracy between different vacua is lifted by other non-perturbative effects [7-10]. 
The domain wall annihilation and the emitted gravitational waves have been extensively 
studied in the literature [11-14]. 

^ The fact that the axion potential can have multiple approximately degenerate vacua has been exploited 
in context of dark energy [3] as well as inflation [4-6]. 
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In this paper we point ont that the annihilation of domain walls also induces the 
baryon asymmetry of the Universe. Suppose that the axion is derivatively coupled to the 
standard model (SM) quarks and/or leptons, 

= ( 3 ) 

where q is a coupling constant. The time derivative of the axion plays a role of the effective 
chemical potential, which spontaneously breaks the CPT symmetry.^ This enables the 
generation of the baryon or lepton asymmetry in thermal plasma if the baryon or lepton 
number is broken, and this is the so-called spontaneous baryogenesis scenario [16-18]. 
The current to which the axion is coupled does not have to coincide exactly with the 
baryon or lepton current; for instance, it could be a U(l) hypercharge current [18] or a 
Peccei-Quinn current [19]. Such derivative couplings to the baryon and lepton currents 
can also be induced if the axion has an anomalous coupling to the SU(2) gauge fields [17]. 
In this case the chemical potential is induced by sphalerons [20, 21], because a non-zero 
time derivative of the axion generates energy difference between the states with different 
winding number and B + L number. Therefore, the chemical potential is expected to be 
suppressed at T > 10^^ GeV where sphalerons decouple from the cosmic expansion. Note 
that there is no such suppression of the effective chemical potential if one starts with the 
derivative couplings with baryon and/or lepton current (more precisely, B — L current), 
as we shall do below. We shall see that, if the axion has such derivative couplings, a 
sizable baryon asymmetry can be generated when the axion domain walls annihilate. 

Before going into details, let us give a rough sketch of our scenario. For simplicity, 
we assume that only two vacua, oi and 02 with ai < 02 , are populated during inflation, 
leading to formation of domain walls separating the two vacua. Generalization to the 
case of multiple vacua is straightforward. After formation, domain walls randomly move 
around at relativistic speed, collide and annihilate continuously, so that the domain wall 
network show the dynamical scaling behavior [22-25]. Every time a domain wall goes 
through some point in space, the held value of the axion changes either from oi to 02 

^ See Ref. [15] for leptogenesis using explicit (non-dynamical) CPT-breaking interactions. 
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or from 02 to aj. Such transition induces a temporal and local chemical potential for 
baryons or leptons. No net baryon asymmetry is generated by the domain wall dynamics 
in the scaling regime, however, because both transitions occur with an equal probability 
and there is no preference of baryons over anti-baryons. The asymmetry between the 
two vacua becomes important when the domain walls annihilate because of the energy 
bias. Suppose that one of the vacua is energetically preferred, e.g., V{ai) < V(a 2 ). When 
domain walls annihilate, the value of a then decreases from 02 to Oi in a region of the false 
vacuum, which gives a preference to baryons over anti-baryons for a certain choice of the 
couplings. Thus, the axion domain wall annihilation can generate the baryon asymmetry 
of the Universe. 

Our scenario has several advantages. First, it is known that the spontaneous baryoge- 
nesis in the slow-roll regime generically leads to baryonic isocurvature perturbations [26], 
which makes the scenario incompatible with high-scale inflation.^ In our scenario, how¬ 
ever, the baryonic isocurvature perturbations can be signihcantly suppressed, because of 
the scaling property of the domain wall network. In particular, our scenario is consistent 
with large-held inhation, and therefore, the required high reheating temperature can be 
realized more easily. Secondly, the axion held value is kept large inside domain walls, 
which enables a large ehective chemical potential even when the axion mass m becomes 
larger than the Hubble parameter. Without domain walls, the spontaneous baryogene- 
sis would become inefficient when the axion starts to oscillate about the minimum [28]. 
Therefore, the axion domain wall baryogenesis scenario works for a wide range of the 
axion mass and the inhation scale. 

Lastly let us comment on diherences of our scenario from other works. In the thick- 
wall regime of the electroweak baryogenesis, the passage of an expanding bubble wall 
generates a non-zero chemical potential, which leaves net baryon asymmetry in thermal 
plasma based on the spontaneous baryogenesis [17, 18] (see also Ref. [29]). The bubble 

^ It is possible to give the axion a mass of order the Hubble parameter in the spontaneous baryogenesis 
using a flat direction [19], thus avoiding the isocurvature constraint. Also, no isocurvature perturbation 
is induced in the gravitational baryogenesis [27]. 
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walls play a similar role to that of domain walls in our scenario. The difference is that 


the electroweak spontaneous baryogenesis relies on the hrst order phase transition of 


two (or more) Higgs helds, and the sphaleron process is exponentially suppressed in the 
symmetry breaking vacuum. As a result, the estimate of the hnal baryon asymmetry 
requires a precise determination of the critical held value as well as detailed analysis 
of the diffusion process during the phase transition [30]. In our scenario, on the other 
hand, the baryon (or lepton) number violation is operative equally in the two minima. 
Also it relies on the domain wall dynamics of a single axion held, whose behavior is well 
studied with numerical simulations. This makes our scenario relatively simple and robust. 
Recently, the authors of Ref. [31] proposed a scenario where the axion has only anomalous 
coupling to SU(2)i gauge helds. They studied a spatially homogeneous axion held in the 
slow-roll regime, and explored the parameter space of the axion mass and decay constant 
preferred by the string axions. The parameter ranges have an overlap with our scenario. 
One diherence is that we start with derivative couplings of the axion with baryon and/or 
lepton currents. Another is that our scenario relies on the domain wall dynamics, while 
Ref. [31] focused on the homogeneous axion held. 

The rest of this paper is organized as follows. In Sec. II, we briehy review the evolution 
of axion domain walls. We estimate the baryon asymmetry induced by the domain wall 
annihilation in Sec. III. The last section is devoted to discussion and conclusions. 

II. AXION DOMAIN WALLS 

Let us consider an axion whose potential is given by 



( 4 ) 


where m and / are the mass and the decay constant of the axion a. We assume that 
two adjacent minima, oi = 0 and 02 = 27r/, are populated with more or less equal 
probability during inhation, and that domain walls separating the two minima are formed 
after inhation. This is the case if the quantum huctuations of the axion, 5a ~ Hini/2Tr, 
is comparable to the decay constant, or if the initial position of the axion is sufficiently 
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close to the local maximum, Omax = tt/. Our scenario can be straightforwardly applied to 
the case in which more than two minima are populated. 

The domain wall solution in a flat spacetime is given by 

ttdwit, x) = 4/ tan“^ exp [m'y{x — vt )~\, (5) 

where x is the spatial coordinate perpendicular to the domain wall, v is the domain wall 
velocity and 7 is the relativistic factor dehned by 7 = 1/a/ 1 — The above solution is 
valid if the thickness of the domain wall ~ m~^ is much smaller than the Hubble horizon, 
i.e., H, where H is the Hubble parameter. The energy density of the domain wall 

is characterized by the tension cr, 

a = 8m/^, (6) 

for the potential (4). 

The domain walls are formed when H m. According to the numerical and analytic 
calculations [22-25], within a few Hubble time after the formation, the domain walls 
quickly follow the scaling law, i.e., 

Pdw ~ crH, (7) 

where there are only one or a few domain walls in each Hubble horizon. The domain 
walls must annihilate and disappear before they start to dominate the Universe, since 
otherwise the Universe would be too inhomogeneous. We assume that there is an¬ 
other shift-symmetry breaking term which generates a bias between the two minima, 
e = V{a 2 ) — V{ai). Then domain walls annihilate rapidly when the energy density of 
domain walls becomes comparable to the energy bias [ 8 - 10 ], 


Pdw ^ 


( 8 ) 


Marginally relativistic axion particles with a typical momentum, k m, are copiously 
produced through the axion domain wall annihilation. Those axion particles soon become 
non-relativistic due to the cosmic expansion [12-14]. In addition, axion coherent oscilla¬ 
tions are produced at the domain wall formation, and we shall discuss their cosmological 
impact later in this paper. 
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The axion particles eventually decay into SM particles through their couplings with 
the SM sector. In general, the axion can have derivative couplings to fermions like (3), 
which are allowed by the shift symmetry (1). Specifically we focus on the case in which 
the axion has derivative couplings only to the SM left-handed lepton currents,"^ 


£ 3 ¥ E (9) 

Our results remain practically unchanged even if one adds additional derivative couplings 
to other SM fermions. If the axion is coupled to the SM sector only through the above 
interaction (9), it mainly decays into a pair of SU(2)i gauge bosons and hypercharge 
gauge bosons through its anomalous couplings [16] . The decay width into a pair of gauge 
bosons is approximately given by 


T 


a 




a'2 \ Njm^ 
2567r3 ^ 10247r3 ) P 


( 10 ) 


where a 2 and a' are respectively SU(2)i and U(l)y gauge coupling constants and Nf is 
the number of of generation, and we will set Nf = ?>m the following. Approximating that 
this is the main decay channel, the axion decay temperature is 


Ta^Sx 10^ GeV 


m y/Yl0^5GeV\ 

IQiiGevJ V f ) 


( 11 ) 


where we have defined the decay temperature by 3H{Ta) = Tq. If those axion particles 
dominate the Universe before the decay, there will be an extra entropy production by the 
axion decay, which dilutes pre-existing baryon asymmetry by some amount. As we shall 
see, the entropy dilution becomes important for a large decay constant and a small axion 


mass. 


^ In a supersymmetric theory, this type of coupling arises from the Kahler potential K — j{A + A^)L^L, 
where A and L are respectively the axion and the lepton supermultiplet, and the lowest component of 
A is given by the saxion and axion as A = s -t ia. 
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III. BARYOGENESIS BY DOMAIN WALL ANNIHILATION 


A. Analytical estimate of the asymmetry 

Now let us discuss baryogenesis by the axion domain walls under the existence of 
the derivative coupling to the lepton current given by (9).® As previously noted, if d is 
non-vanishing, the derivative couplings behave like an effective chemical potential, 

^ 3 “ = fled” + . . . . (12) 

where fies = d/f is the effective chemical potential for the lepton number (L). 

The axion domain walls can generate the effective chemical potential because of the 
large spatial gradient of the axion held inside the wall. Since domain walls are moving at 
nearly the speed of light, the time derivative of the axion held at some hxed spatial point 
becomes large while domain walls are passing through. The ehect of the gradient term is 
negligible if the domain wall is sufficiently thick compared to the dihusion length. 

If the L-number violating operator is in equilibrium, and if the chemical potential is 
spatially homogeneous, the diherence of number densities between lepton and anti-leptons 
would be produced as ~ for /leff T, where we have taken into account 

the spin degrees of freedom and the number of generation. It depends on the rate of 
the L-number violating process as well as the domain wall dynamics whether the lepton 
asymmetry reaches the equilibrium value in the expanding Universe. One needs to solve 
the Boltzmann equation for the lepton asymmetry, hl = n£ — ni, 

rii + SHul = -T{nL - (13) 

where T is the interaction rate for the L-violating processes. Note here that the chemical 
potential in depends on the position and velocity of domain walls. 

As the L-number violating operator, we consider AL = 2 scattering processes, i£ y-)- 
HH, iH yy iH, which are mediated by heavy right-handed Majorana neutrinos in the 

® Instead, one may use an anomalous coupling of the axion to the SU(2) gauge fields, in which the baryo¬ 
genesis works similarly as long as the sphalerons are in equilibrium at the domain wall annihilation. 
(See the discussion in Sec. I). 




seesaw mechanism [32-35]. Here and in what follows we assume that the right-handed 
neutrinos are so heavy that they can be integrated out in our analysis. The interaction 
rate for the AL = 2 processes is roughly given by [36] 


I rsj - 

'jrS 

^ ^EW 

where vew = 174 GeV and rrii with i = 1, 2, 3 denotes the mass of three active neutri¬ 
nos. The decoupling temperature of the L-violating process in the radiation dominated 
Universe is 



Tdec ~ 3 X 10^^ GeV, 


(15) 


where we have assumed the normal ordering for the neutrino mass differences and used 
the experimental value, ^ ~ ~ 2.4 x 10“^ eV^. For the reheating temperature 

Tr lower than Tdec, the T-violating process remains decoupled from the cosmic expansion. 
As we shall see below, even in this case, a non-zero lepton asymmetry is induced by the 
domain wall annihilation. 

Let us hrst consider an ideal situation where a domain wall passes through the origin 
X = 0 at t = t^w- Using Eq. (5), the effective chemical potential at the origin evolves 
with time as 

cosh[m 7 u(t — tow)] 

It takes roughly At ~ {m'yv)~^ for the domain wall to pass through the origin, and so, 
the induced lepton asymmetry by passage of the domain wall is estimated as 


riL ~ Tn^At ~ TTl 


(17) 


Note that the lepton asymmetry becomes independent of the velocity of the domain walls. 
As the domain wall passes through, a similar amount of the lepton number density will 
be induced inside the Hubble horizon. 

In the scaling regime, domain walls randomly move around inside the Hubble horizon 
so as to collide and annihilate continuously. In particular, since there is no preference 
for either of the vacua, the effective chemical potential can be positive or negative with 
equal probability. Therefore there will be no net lepton asymmetry left, even though 
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some amount of the lepton asymmetry with either positive or negative sign is induced 
each time a domain wall passes through. Such lepton asymmetry has fluctuations of order 
unity inside the Hubble horizon, but it has no sizable fluctuations at superhorizon scales, 
because of the scaling property of the domain-wall network. 

A non-zero net lepton asymmetry is induced when domain walls annihilate and dis¬ 
appear owing to the energy bias. This is because one of the two vacua is energetically 
preferred, inducing an effective chemical potential with a hxed sign in the false vacuum 
which occupies about half of the space. Again, the scaling property of the domain wall 
network ensures that there is no isocurvature perturbations at super-horizon scales. 

The hnal lepton asymmetry is generated within about one or a few Hubble time before 
the domain wall annihilation. In particular, the maximal possible value of the lepton 
asymmetry is obtained when the domain wall annihilation takes place at the decoupling 
of the L-violating processes. The reason is as follows. If the domain wall annihilation takes 
place before the decoupling of the L-violating processes, the lepton asymmetry induced by 
the domain wall annihilation will be washed out. On the other hand, if the domain wall 
annihilation occurs after the decoupling, the induced asymmetry tends to be suppressed 
because the L-violating process is inefficient. The maximum asymmetry is therefore 


s 




max 


45 T 
T^‘^g*s T 


-10 


-6 


dec 


( 18 ) 


where s and g^s are respectively the entropy density and the relativistic degrees of freedom. 
We have substituted g^,s = 106.75 and the decoupling temperature (15) in the second 
equality, assuming the radiation-dominated Universe. The negative sign is inserted in the 
second equality to obtain positive baryon asymmetry through sphalerons. 

If the reheating temperature Tr is lower than the decoupling temperature Tdec, the 
interaction rate for the L-violating processes never exceeds the expansion rate of the 
Universe. One can see this by noting that T/H reaches the maximal value (smaller 
than unity) at the reheating as long as the temperature of the dilution plasma obeys 
T ~ before the reheating. Hence the maximal asymmetry in this case is 

obtained if the domain wall annihilation occurs at the reheating, and it is roughly given 
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by 


2 


s 


R 


s 


T, 


R 


z 


dec 


(19) 


We shall see later in this section that the maximal asymmetry is indeed generated if the 
domain wall annihilation takes place at T = min[Tdec, 7 r]- 


B. Necessary conditions for successful baryogenesis 

Here let us discuss some necessary conditions for the successful domain wall baryo¬ 
genesis. First, the domain wall dynamics should have negligible back reaction from the 
generated lepton asymmetry in the plasma. As the domain walls move in the plasma, 
some amount of the lepton asymmetry is induced because of the effective chemical po¬ 
tential (16). The interaction with the generated asymmetry induces a back reaction, 
which would act as a frictional force on the domain wall dynamics. The back reaction is 
negligible, and the domain walls follow the scaling law if 

aH > /iesriL (20) 

at the domain wall formation (ifform ~ where is given by (17). 

Secondly, the domain wall must be sufficiently thick to justify our analysis where we 
have neglected dissipation of the asymmetry. The thickness of the wall is roughly m~^ and 
the typical mean free path of the particle in plasma is of order T~^. Thus, the thick-wall 
condition is given by 

Tlnn > rn, (21) 

where Tann denotes the temperature at the domain wall annihilation. 

Thirdly, we have assumed that the domain wall annihilation takes place well after the 
domain wall network start to follow the scaling law. It takes a few Hubble time after the 
formation to reach the scaling regime, and therefore we conservatively require 

Hform ~ m, > lOTTann) (2^) 

where inform and Ffann are the Hubble parameter at the domain wall formation and anni¬ 
hilation, respectively. 
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Fourthly, we require that the decay constant is larger than the quantum fluctuations 
of the axion to ensure the validity of analysis using the potential (4). Specifically, we 
impose a lower bound on / as 


/ > 


H; 


inf 


271 ’ 


(23) 


where iifinf is the Hubble parameter during inflation. If this bound is not satisfied, the 
corresponding U(l) symmetry may be restored, or the saxion held may be destabilized. 

Finally we assume that there is (effectively) only single path connecting the two vacua 
tti and 02 . Apparently this is not satisfied if a U(l) symmetry is explicitly broken down to 
Z 2 . In this case there are two paths (clockwise and counter-clockwise) connecting the two 
vacua. In other words, there appear two kinds of domain walls with the same number. 
This can be understood by noting that the two types of the domain walls are attached to 
cosmic strings associated with the spontaneous break down of the U(l) symmetry. If the 
tensions of the two type of domain walls are equal, they would start to annihilate at the 
same time and sweep equal spatial volume with positive and negative chemical potential, 
resulting in no net baryon asymmetry. On the other hand, if there is an explicit breaking 
of the Z 2 symmetry such that one type of domain walls has a larger tension than the other 
one, the domain walls with a smaller tension would start to annihilate first by the energy 
bias between the two vacua and sweep a larger spatial region, producing a net baryon 
asymmetry. Therefore, our scenario works even if there are multiple paths connecting 
the two vacua (namely if there are multiple types of domain walls), as long as one of the 
multiple paths is energetically favored. If there are multiple vacua, or if the symmetry is 
non-linearly realized, our scenario works by a similar argument. 

In the numerical calculations we impose the above conditions to ensure successful 
domain wall baryogenesis. It turns out that all the conditions are easily satisfied for the 
parameters of our interest. 
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C. Numerical calculations 


The net lepton asymmetry is effectively induced by the domain wall annihilation, 
during which domain walls sweep typically about a half of the space. To model the domain 
wall dynamics during the annihilation, we approximated the situation by a single domain 
wall passing through the origin, where we numerically solve the Boltzmann equation (13), 
combined with the evolution equations for the energy density of the inflaton (p/) and 
radiation (p^), 

Pi+ ?>Hpi =-Tipi, pr + 4:Hpr = Tipi, (24) 

where T/ is the decay rate of the inflation, and we define the reheating temperature in 
our analysis by 3H{Tii) = Tj. This approximation is valid because no net asymmetry 
is induced during the scaling regime, and so, we can focus on the domain wall dynamics 
during the one or a few Hubble time before the annihilation. 

In Fig. 1 we show the induced lepton asymmetry as a function of the domain wall 
annihilation temperature for various values of the reheating temperature. In the top and 
bottom panels, we have set the axion mass to be m = 10^^ GeV and 10^^ GeV, respectively. 
Here we have not taken into account the entropy production by the subsequent axion 
decay, which we shall return to in a moment. As expected, the maximal asymmetry is 
obtained when Tann — min(Tdec, Tr), in good agreement with the analytic estimate (18). 
In the bottom panel, one can see that the lepton asymmetry is highly suppressed in the 
case of e.g. Tann > Tdec and Tr = 10^“^ GeV. This is because the asymmetry induced 
by the domain wall annihilation is subsequently washed out by the T-number violating 
processes in equilibrium. In general, we expect that the wash-out process is efficient when 
Tr > Tann > ^dec- 

At the domain wall annihilation, marginally relativistic axions are copiously produced, 
and they may come to dominate the Universe before they decay into gauge bosons. Once 
the axion dominates the Universe, its subsequent decay produces a large entropy, diluting 
pre-existing asymmetry. Thus, the hnal baryon asymmetry is hxed after the axion decay. 
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(a) m = 10^^ GeV 



(b) m = 1012 GeV 

FIG. 1: The induced lepton asymmetry as a function of the domain wall annihilation temper¬ 
ature for various values of Tr and the axion mass m = 10^^ GeV (top) and 10^^ GeV (bottom). 
The vertical dotted (magenta) line represents the decoupling temperature of the L-number 
violating processes in a radiation-dominated Universe. Note that the subsequent entropy dilu¬ 
tion by the axion decay is not taken into account here. We have imposed the condition (22), 
m > lOFfami) which corresponds to the left end point of each curve. 
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if there is entropy dilution. Taking into account the sphaleron process®, the resultant 
baryon asymmetry is estimated as 


Ub 

s 



1 . 

- X A X 
2 


s 


(25) 


where A is the dilution factor by the axion decay given by 

(DW annihilation before reheating) 

(DW annihilation after reheating) 


A = 


min I 


mini 1, 

» (jHann 


(26) 


The numerical factor 1/2 comes from the fact that the transition from the false vacuum 
to the true vacuum takes place in about half of the whole space. 

In Figs. 2 and 3 we show the contours of the hnal baryon asymmetry, n^/s, in the m-f 
plane for various values of Tr. Here we have set = min(Tdec,TR) so that the baryon 
asymmetry takes the largest possible value for a given reheating temperature. The baryon 
asymmetry can be suppressed by either increasing or decreasing Tann (see Fig. 1(a)). One 
can see that a sufficient amount of baryon asymmetry, ub/s > 10“^®, can be generated 
for Tr > 2 X 10^^ GeV. In the lower shaded (magenta) region, there is no entropy dilution, 
i.e., A ~ 1, and so, hb/s takes a constant value. As / becomes large, hb/s decreases 
owing to the entropy dilution factor A ^ 1. This is because, as / increases, the energy 
density of the axion particles increases and the lifetime of the axions becomes longer. The 
horizontal dashed (green) lines and dash-dotted (cyan) lines represent the lower bound 
on the axion decay constant, / > ifinf/27r, for = 10^"^ GeV and aH > /refr^^L|H=m, 
respectively (cf. (20) and (23)). The yellow-shaded region in upper right corner in Fig. 3 
is ruled out from the domain wall domination at annihilation. Below the dotted (blue) 
line, baryonic isocurvature perturbations and their non-Gaussianity would exceed the 


® We have neglected the sphaleron effects during the domain wall annihilation, for simplicity. This 
approximation is valid for most of the parameters of our interest, because, as we shall see, successful 
baryogenesis requires Tann ^ 2 x 10^^ GeV, while the sphalerons are decoupled at T > 10^^ GeV. Even if 
sphalerons are in equilibrium at the domain wall annihilation, the resultant baryon asymmetry changes 
only by a factor of 0(1), and our main results remain valid. 
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observational bound, if the L-number violating rate (14) is valid at the domain wall 
formation. In other words, in the region slightly below the dotted (blue) line, baryonic 
isocurvature perturbations and their non-Gaussianity may be found in the near future 
observations. We will discuss this issue in the next subsection. 

D. Baryonic isocurvature perturbations 

Here we discuss baryonic isocurvature perturbations in our scenario. Here we do not 
distinguish lepton asymmetry and baryon asymmetry, as we are concerned with the hnal 
baryon asymmetry at the CMB epoch. First, let us consider baryon asymmetry generated 
by the domain wall annihilation, FDw,ami = '^s/s|DW,ann- As domain walls are spatially 
localized objects, lDw,ann has initially large spatial fluctuations of order unity at subhori¬ 
zon scales. Such small-scale fluctuations asymptote to zero in the course of evolution, 
because of diffusion processes of quarks and leptons. At super-horizon scales (e.g. the 
CMB scales), on the other hand, how.ann has no isocurvature fluctuations because of the 
scaling property of the domain wall network. This results stand in sharp contrast to the 
usual spontaneous baryogenesis in the slow-roll regime [26]. 

Secondly, we turn to baryon asymmetry generated right after the domain wall forma¬ 
tion. We have assumed that the axion acquires sufficiently large quantum fluctuations 
during inflation so that the two adjacent vacua are realized randomly in each Hubble 
horizon. This leads to the formation of domain walls when the Hubble parameter be¬ 
comes comparable to the axion mass, if ~ m. At the same time, the axion coherent 
oscillations are induced. The dynamics of axion coherent oscillations, especially its mo¬ 
tion in the slow-roll regime, generates the baryon asymmetry in the background thermal 
plasma as in the usual spontaneous baryogenesis. Let us denote the baryon asymmetry 
by Hose- As the axion has initially large quantum fluctuations at super-horizon scales, Wsc 
has isocurvature fluctuations at large-scales, which is the counter part of the baryonic 
isocurvature fluctuations in the spontaneous baryogenesis in the slow-roll regime. In our 
case, the size of the baryonic isocurvature perturbations, ^Wsc/Wsc, is expected to be of 
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(a) Tr = 2 X 10“ GeV 



(b) Tr = 10“ GeV 


FIG. 2: Contours of the final (maximal) baryon number asymmetry in the m-f plane for 

Tr = 2x 10^^ GeV (top panel) and 10^^ GeV (bottom panel). We assume Tann = min[rR, Tjec] so 

that the baryon asymmetry becomes maximal. The solid (red) lines correspond to the contours 

of ur/s = from top to bottom. In the shaded (magenta) regions, there is no 

entropy dilution (i.e. A = 1), and ur/s takes a constant value ur/s = 8.5 x 10“^^ (top panel) 

and 2.1 X 10“® (bottom panel). In the cyan-shaded region, the thick wall condition is violated. 

Baryonic isocurvature perturbations and their non-Gaussianity will be too large below the dotted 

(blue) line, as long as one extrapolates the L-violating interactions to the domain wall formation. 

See the text for discussion on this issue. The horizontal dashed (green) lines represent the lower 

bound on /, / > do ~ FIinf/27r for = 10^^ GeV. 
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(a) Tr = 1013 GeV 



FIG. 3: Same as Fig. 2 but for Tr = GeV (top panel), GeV (bottom panel). The 
solid red lines correspond to the contours of ub/s = from top to bottom and the 

magenta shaded regions correspond to the maximal value, 1.6 x 10“^ (top panel) and 6.9 x 10“^ 
(bottom panel). The dash-dotted cyan line represents the lower bound from the back reaction 
and the yellow-shaded region is ruled out from domain wall domination. 


order unity. This can be understood by noting that the chemical potential can be either 
positive or negative, depending on which vacuum the axion is rolling down to. After the 
commencement of oscillations, the scalar wave dynamics between walls are random and 
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complicated. In particular, the spatially averaged effective chemical potential is zero, and 
no fluctuations at super-horizon scales are induced by the dynamics in the scaling regime. 
Therefore, l^sc and its fluctuations at large scales receive the main contribution from the 
domain wall formation when H ^ m. 

Finally, the domain-wall dynamics toward the scaling regime will also induce the baryon 
isocurvature perturbations. For domain walls to be formed, or more precisely, for inhnitely 
long domain walls to be formed, the probabilities to realize the two vacua must be com¬ 
parable, but they do not have to be exactly equal to each other. It implies that, when 
domain walls are formed, the spatial volume of one of the vacua is generically larger (or 
smaller) than that of the other by (at most) a few tens of percent. The ratio of the two 
volumes will quickly converge to unity as the domain-wall network approaches the scaling 
evolution. This is because the two vacua are degenerate in energy and there is no prefer¬ 
ence to one over the other once the scaling regime is reached. In this process toward the 
scaling regime, there is an overall transition from one of the vacua to the other, which sim¬ 
ilarly induces the baryon asymmetry. Let us denote the asymmetry by low, form- As the 
bias of the spatial volumes is induced by the quantum fluctuations of the axion, low, form 
has isocurvature fluctuations at large scales. The magnitude of low, form is expected to 
be comparable to Fqsc, and the sign is opposite. So, there is a partial cancellation, but in 
general, there is no exact cancellation. For our scenario to work, both Fosc and low, form 
must be sufficiently suppressed, since otherwise the baryonic isocurvature perturbations 
and their non-Gaussianity, would be too large to be consistent with observations. 

The baryon asymmetry generated at the domain-wall formation can be suppressed as 
follows. If the lepton-number violation processes are in equilibrium between the formation 
and annihilation of domain walls, the initial asymmetry Tosc and low, form can be washed 
out. This is the case if the reheating temperature is higher than 10^^ GeV. For lower 
reheating temperature, the lepton-number violating processes remain decoupled all the 
time. Then, IGsc and low, form can be suppressed if the lepton-number violating rate is 
much smaller than the Hubble parameter at the domain formation. 

In our numerical calculations, we have estimated |T()sc| ~ I how, form I by following the 
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motion of a test domain wall which goes through a fixed position aX H = m. Using the test 
domain wall as background classical held evolution, we have calculated the induced baryon 
asymmetry in the plasma by solving the Boltzmann equation. By doing so, we effectively 
evaluate |low,form| (or |Uosc|) at the formation, neglecting the complicated dynamics of 
the scalar waves and domain-wall evolution, which do not have any preference to baryons 
over anti-baryons. 

The current constraint on the matter isocurvature perturbation S from the Planck 
observation reads Vs < 8.7 x 10“^^ [37]. Using the fact that baryon isocurvature per- 
turbation is written as ~ Sflb/^m — 0.15((5r2b/r2b), we obtain the constraint on the 
baryon isocurvature perturbations as 6Qb/^b ^ 6 x 10“^. Since the baryons produced 
by the axion coherent oscillations or domain wall dynamics toward the scaling regime is 
0(1) in the present scenario, flb,osc/^b ^ 6 x 10“^ must be satished in order to avoid too 
large isocurvature perturbations. Then, we obtain the constraint on the resultant baryon 
asymmetry induced by the coherent oscillations, 

< 5 X 10-15, (27) 

s s ilb 

and a similar bound on the asymmetry induced by the domain wall dynamics toward the 
scaling regime. This upper bound is shown by a dotted (blue) line in Figs. 2 and 3. 

The baryon isocurvature perturbations may be further suppressed in some particular 
situations. For example, one can consider a case in which the U(1 )b_l gauge symmetry 
is still unbroken the onset of the axion oscillation and it gets spontaneously broken before 
the domain wall annihilation. In such a case, there is no lepton number violating operators 
and no baryon asymmetry is induced until the spontaneous break down of the U(1 )s_l 
symmetry. If the domain wall network already follows the scaling law when the U(1 )s_l 
symmetry gets spontaneously broken, no baryon isocurvature perturbation is generated 
by the coherent oscillations or domain wall dynamics. Interestingly, cosmic strings are 
formed after the spontaneous breaking of U(l) 5 _i;, and they can emit a sizable amount of 
gravitational waves which can be within the reach of future observations [38]. 
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IV. DISCUSSION AND CONCLUSIONS 


Collapsing domain walls are cosmological sources of gravitational waves [11-13]. The 
gravitational wave spectrum is peaked at a frequency, 

/p.., -1« kH. . (28) 

corresponding to the Hubble horizon scale at the domain wall annihilation [14]. Here ^ 
and Tx are defined as 

^ = min ^1, T^nn)- (29) 

For successful baryogenesis, Tx must be higher than 2 x 10^^ GeV, and so, the peak 
frequency is at (9(100) kHz or higher, which is too high to be detected by near future 
observations. We note however that there have been proposed several new detection 
techniques with the sensitive frequency region around MHz [39, 40], which may be able 
to probe gravitational waves produced in our scenario. 

So far, we have considered the L-number violating processes mediated by heavy right- 
handed neutrinos in the seesaw mechanism. Other types of the baryon/lepton vio¬ 
lating operator is also possible and the corresponding decoupling temperature for the 
baryon/lepton violating processes could be lowered. One of the examples is the R-parity 
violating operator, 

^ ~ 2^ijkLiLjEk (30) 

in the supersymmetric Standard Model. In this case, the interaction rate for the L- 
violating processes scales as F cx T® for T ^ and F cx T for T S> where is 
the slepton mass. For instance, if we take A ~ 10“® and > 10® GeV, the L-violating 
process marginally reaches equilibrium and soon decouples at T^ec 10® GeV. Since the 
maximal possible value of lepton asymmetry is roughly given by ul/s ~ O-lT^ec/Mp from 
the first equality in (18), successful baryogenesis is possible with Tann ~ 10® GeV. In this 
case, the peak frequency of the gravitational waves from the domain wall annihilation can 
be within the sensitivity range of the ground-based detector such as advanced-LIGO [41] 
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and KAGRA [42, 43]. For instance, if we take Tji ^ m ^ 10® GeV and / ~ 10^^ GeV, 
domain walls dominate the Universe at the annihilation and the peak freqnency falls in 
the sensitivity range of these experiments. A naive order-of-magnitnde estimate snggests, 
however, that the signal strength is a few orders of magnitnde smaller than the predicted 
sensitivity, and either some deviation from the scaling regime or fnrther improvement of 
the sensitivity wonld be necessary to directly probe snch signals. 

In this paper we have proposed a baryogenesis scenario nsing axion domain walls. Ax- 
ion domain walls are prodnced if the axion acqnires snfficiently large qnantnm flnctnations 
during inflation or if it initially stays sufficiently close to the local maximum. While no 
net baryon asymmetry is produced in the scaling regime, collapsing axion domain walls 
produce a large enough baryon asymmetry to explain the observed value. This is because 
the energy bias between the two vacua, and therefore between baryons and anti-baryons, 
becomes relevant only when domain walls annihilate. In particular, baryon isocurvature 
perturbations can be signihcantly suppressed in our scenario, either because the asymme¬ 
try produced by the initial field conhgurations is washed out by the L-number violating 
interactions in equilibrium, or because the L-number violating interaction is simply sup¬ 
pressed at the domain wall formation. In some parameter region, baryon isocurvature 
perturbations and their non-Gaussianity are suppressed, but non-negligible, which may 
be detected by future observations. Our scenario works together with high-scale inflation 
which predicts a large tensor-to-scalar ratio within the reach of future B-mode observa¬ 
tions. The required relatively high reheating temperature can be realized in high-scale 
inflation more easily. This should be contrasted to other spontaneous baryogenesis scenar¬ 
ios in which the inflation scale is severely constrained by the isocurvature perturbations. 
Although we have focused on the axion domain wall throughout this paper, our analysis 
can also be straightforwardly applied to a wide class of domain walls such as the Standard 
Model Higgs domain wall [44, 45]. 
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